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Abstract 

In this work, we construct the non-relativistic Lee model on some class of three dimen- 
sional Riemannian manifolds by following a novel approach introduced by S. G. Rajeev 
PP| . This approach together with the help of heat kernel allows us to perform the renor- 
malization non-perturbatively and explicitly. For completeness, we show that the ground 
state energy is bounded from below for different classes of manifolds, using the upper 
bound estimates on the heat kernel. Finally, we apply a kind of mean field approximation 
to the model for compact and non-compact manifolds separately and discover that the 
ground state energy grows linearly with the number of bosons n. 

1 Introduction 

The Lee model, originally introduced in [2], is an exactly soluble and renormalizable 
model that includes the interaction between a relativistic bosonic field ("pion") and a 
heavy source ("nucleon") with one internal degree of freedom, which has two eigenvalues 
distinguishing "proton" and "neutron" . By heavy, we mean that the recoil of the source 
is neglected. Although this model is not very realistic, it reflects important features of 
nucleon-pion system and presents a powerful aspect that one can do the renormalization 
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without the use of perturbation techniques. Moreover, the complete non-relativistic ver- 
sion of this model that describes one heavy particle sitting at some fixed point interacting 
with a field of non-relativistic bosons is as important as its relativistic counterpart. It is 
simpler than its relativistic version because it is possible to renormalize the Hamiltonian 
of the system with only an additive renormalization of the mass (energy) difference of the 
fermions. It has been studied in a textbook by Henley and Thirring for small number of 
bosons from the point of view of scattering matrix [3] and there are further attempts in 
the literature from several point of views [4j. It is possible to look at the same problem 
from the point of view of the resolvent of the Hamiltonian in a Fock space formalism 
with arbitrary number of bosons (in fact there is a conserved quantity which allows us to 
restrict the problem to the direct sum of n and n + 1 boson sectors). This is achieved in a 
very interesting unpublished paper by S. G. Rajeev [TJ, in which a new non-perturbative 
formulation of renormalization has been proposed. We are not going to review the ideas 
developed in there. Instead, we suggest the reader to read through the paper [I] to make 
the reading of this paper easier. 

Following the original ideas developed in [I], we wish to extend the non-relativistic 
Lee model onto the Riemannian manifolds with the help of heat kernel techniques, hoping 
that one may understand the nature of renormalization on general curved spaces better. 
In this work, for the sake of simplicity we ignore the motion of the heavy particle and 
take its position a as a given fixed point on the manifold. The construction of the model 
is simply based on finding the resolvent of the regularized Hamiltonian H t and show that 
a well-definite finite limit of the resolvent exists as e — > + (called renormalization) with 
the help of heat kernel. We prove that the ground state energy for a fixed number of 
bosons is bounded from below, using the lower bound estimates of heat kernel for some 
class of Riemannian manifolds, e.g., Cartan-Hadamard manifolds (also explicitly H 3 ), 
the minimal submanifolds of M 3 and closed compact manifolds with nonnegative Ricci 
curvature. We also study the model in the mean field approximation for compact and 
non-compact manifolds separately and prove that the ground state energy grows linearly 
with the number of bosons for both classes of manifolds. 

The paper is organized as follows. In the first part, we construct the model and show 
that the renormalization can be accomplished on Riemannian manifolds. Then, we prove 
that there exists a lower bound on the ground state energy. Finally, the model is examined 
in the mean field approximation. 



2 The Construction of the Model 



We start with the regularized Hamiltonian of the non-relativistic Lee model on a three 
dimensional Riemannian manifold (Ai,g) with a cut-off e. Adopting the natural units 
(h = c = 1), one can write down the regularized Hamiltonian on the local coordinates 
x = (xi,x 2 ,x 3 ) G M 

H e = H + H l!€ , (1) 



where 



H = Jd g x<l>\x)(-±V* + m)<l>(x), (2) 
H Ij6 = //(e) 1 ° Z + A I d g x p e (x,a) ((f)(x)a- +(fJ(x)a + ^ . (3) 
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Here, d g x = det gij dx is the volume element and V 2 is Laplace-Beltrami operator or 
simply Laplacian, and <^(x), </>(#) is the bosonic creation-annihilation operators defined 
on the manifold with the metric structure gij. Sometimes we shall write (j> g (x) in order 
to specify which metric structure it is associated with but for now we simply write down 
4>(x). Also, p e (x,a) is a family of functions which converge to the Dirac delta function 
5 g (x, a) (with the normalization J M d g x 5 g (x) = 1) around the point a on Ai as we take 
the limit e — ► + . The Pauli spin matrices a± = ^(ci ± ia-i) and 03 are regarded as a 
matrix representation of the fermionic creation and annihilation operators acting on C 2 
and /t(e) is a bare mass difference between the "proton" and "neutron" states of the two 
state system ("nucleon"). Its explicit form will be determined later on. Although the 
number of bosons is not conserved in the model, one can derive from the equations of 
motion that there exists a conserved quantity 

Q = - 1 + / d g x 4>\x)<t)(x) 
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which takes only positive integer values. If Q = n £ Z + , we have spin- up state ("proton") 
with n bosons or spin-down state ("neutron") with n — 1 bosons. We can think of the 
latter as a bound state of the system. If one considers the Hamiltonian without the cut-off 
e, it can be shown that the bound state energy diverges. Before discussing how to deal 
with the infinities, we must see how the infinities emerge in our model. The simplest way 
to realize this is just to look at the sector that contains the "neutron" or the "proton" 
with one boson, which corresponds to Q = 1, that is, we propose the following eigenstate 
ansatz [5]: 

k ^ )=| '/^(^M|o) (4) 

u\0) 



For simplicity we explicitly perform our calculations for compact manifolds here, but 
our result is also valid for non-compact manifolds, which we are interested in. If the 
manifold M is compact, then the Laplacian has a discrete spectrum and there is a family 
of orthonormal complete eigenfunctions f a {x) G L 2 (M.) [6] satisfying, 
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-V 2 g Ux) = af a (x), 

d g X fa( X )fv'(x) = , 

J2fZ(x)f a (y) = S g (x,y). (5) 



In fact, one can extend these expressions to some noncompact manifolds, and they also 
satisfy these properties by an appropriate generalization of the measures to the continuous 
distributions in the sense of [7j. 

From the eigenvalue equation H\u, ip) = E\u, tjj), we find the set of equations in terms 
of the bosonic wave function tp(x) = \\^, a f a (x)tp(a) 

u\f*(a) 

na) = . r , (6) 

u{E-n) = X^Ma^a). (7) 

If we substitute the equation ([6j) into the equation ([7|), we obtain 

r E - (& + "•) 

Expressing this equation in terms of heat kernel K s /2 m (x,x') = (x|e sV 9 //2m |x / ), we get 

^ = E + \ 2 ds K s/2m (a,a)e- s( - m - E K (8) 
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Let us make a short digression on heat kernels. When we want to emphasize the metric 
structure gtj on which the heat kernel is based, we shall use the notation K s /2 m (x,x';g) 
throughout the paper. Although the notion of heat kernel can be defined on any Rie- 
mannian manifold, the explicit formulas only exist for some special class of manifolds, for 
example, Euclidean spaces M. d [8] and hyperbolic spaces M d [9]. Some of the well known 
properties of the heat kernel on A4 in dimensionless parameters t and x are 

Kt(x,y) = Kt(y,x) Symmetry property , 

dKAx,x') , 

\/ a K t {x,x) = Heat equation , 

lim Kt(x,x') = 5 g (x,x') Initial condition , 

/ d g x K tl {x,z)K t2 (z,y) = K tl+t2 (x,y), Reproducing property , (9) 
JM 



for t > only. If in addition Ai is a compact manifold, we have 



K t (x,y) = Y,e- ta f;(x)Uy) 



(10) 
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which converges uniformly on M x A4 . Again we assume that this has a proper analogous 
expression for non-compact manifolds we are interested in. When the manifold A4 is a 
complete Riemannian manifold with Ricci curvature bounded from below then the heat 
kernel satisfies the stochastic completeness property [8j [10] : 



On a compact manifold stochastic completeness is always satisfied [llj . 

The integral in the equation ([8]) diverges due to the asymptotic expansion of the 
diagonal part of heat kernel near s = for any three dimensional geodesically complete 
manifold whose injectivity radius has a positive lower bound [12j 



Here the Uk{a, a) are functions given in terms of curvature tensor of the manifold and its 
covariant derivatives at the point a. As a result of this, the bound state energy becomes 
divergent. In flat spaces, one can do the similar calculations in momentum space and find 
that [5] 



where uj(p) = £— +m. This momentum integral blows up at large values of momentum in 
three dimensions. The problem is basically taking the integral over all momenta because 
our no-recoil approximation breaks down for large enough momenta. So, we introduce an 
ultraviolet cut-off A in the upper bound of the integral. Since large momenta means small 
distances, this cut-off corresponds to putting a small distance cut-off in coordinate space. 
Performing the calculations in coordinate space one see that small distance cut-off can 
be replaced with a short "time" cut-off e in the lower limit of the integral (jSJ. Here we 
show that the idea of short "time" cut-off will work on Riemannian manifolds, whereas 
the momentum cut-off is not a natural method to use. 

Therefore, we first regularize the Hamiltonian, that is, introduce the cut-off e on the 
lower bound of the integral and make the parameters in the Hamiltonian depend on e 
such that all physical quantities are independent of it. So we define 







(12) 




(13) 



where E is traded with fi which is defined as the physical energy of the composite state 
which consists of a boson and the attractive heavy neutron at the center. Now, using ©, 
([7]) and (fT3|) . we get the finite expression for the bound state energy 



E = fi + A / ds K s / 2m (a, a) 
Jo 



e -s(m-^) _ e -s(m-E) 



(14) 



Here it is easy to see that E = \x is a possible solution to this equation. Then, one can 
calculate the bosonic wave function tp(x) for E = fi 

/•oo 

4(x) = -u\ / ds K s/2m (x, a)e- s( - m -^ , (15) 
J o 

if x a. Although the wave function is divergent as x — > a, it is square integrable as we 
will see. If we substitute (|15p into our ansatz @, we get 

V ^ (a) | > 



\uip) = u 



10) 



(16) 



Normalizability of (|16p can be easily seen by using the properties of heat kernel 



A 2 / d g x / d Sl ds 2 K Sl/2m (x,a)K S2/2m (x,a)e-^ +s ^ m -^ 
J M JO 

| I I" dt] K s/2m (a,a) e - s{m -^ = H sds e " s ( m -^ K s/2m (a,a) 
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and result we find the normalization to be 



1 + A z / sds e 
o 



-s(m— n) 



K s /2m(a,a) 



-1/2 



(17) 



(18) 



This integral is finite due to the short and long time behaviour of heat kernel. 

One can also consider the scattering of a boson from the "proton" at rest on a non- 
compact manifold^. The inhomogeneous Schrodinger equation (H — E)\u,ip) = \v,x) 
leads to 



X {a) - \uf*{a) 
^ + m~E 



and 



A]T/»v^) + ^-£) 



(19) 
(20) 



If we substitute CGD into ([20]), we find 

l/»| 2 



A it > — +u(a — E) = v — A > — - 

^-?- + m-E xt ' ^ ^- + m- E 

2m 



2m 



^The same normalization can also be found by writing the operators 0(a) in the eigenbasis f a (a) of the 
Laplacian. 

2 Compact manifolds have only discrete spectrum. 



If we express the above equation in terms of heat kernel, we immediately see that the 
integral is divergent due to singularity near s = 0. So we must do the regularization by 
introducing a cut-off e in the lower limit of the integral and choose //(e) as above. Taking 
the limit e — > + we get the following finite expression 

(poo poo 
X 2 u J ds K s/2m (a, a) e"^"^ - X 2 u J ds K s/2m (a, a) e^™"^ 

= v - X f d g x x(x) ( I ds K s/2m (x, a) e^" 1 "^ 



g x X{X) / us i\ s / 2m \ 
M \J0 



From the above equation, it follows that 
u = u[v,x] 



X 2 I ds K s/2m (a, a) ( e~< m ~^ - e"^^) 



oo 



v - X J ds tj d g x K s/2m (x, a) X (x)j e~< m -^ 
We can also read off ij){x) from the equation ([T9]) when x ^ a 

i>{x) = / ds e-^) / d g yK s/2m (x,y) X (y) - Xu[v, X ] / ds K s/2m (x, a) e~< m -^ , 
Jo J Jo 

(21) 

and ip(a) = A _1 (t> — u[v,x]{n ~ E)) as a result of the equation (f20j) . It is important to 
remind the reader that these expressions should be understood in the sense of analytic 
continuation in the complex S-plane to their largest domain of definition. Indeed if the 
real part of E is smaller than m these integrals all make sense, and the resulting expressions 
are just the Green's function, or the resolvents for the Laplace operator, which exists away 
from the positive real axis. 

Up to now, we have shown that non-relativistic Lee model on a manifold is divergent 
and can be renormalized with the help of heat kernel by considering the problem for Q = 1 
sector. However, it is not clear in this formulation, how we can extend this method to the 
case that have arbitrary number of particles, say Q = n sectors and show that the ground 
state energy is bounded from below. Nevertheless, we have an alternative and powerful 
method which is developed by Rajeev [T]. From now on, we will follow his approach in 
order to construct and develop the model on a manifold for any sector. 

Let us first express the regularized Hamiltonian as a 2 x 2 block split according to C 2 : 

H E= ( Hq-E Xj M d g x p e (x,a)ft(x) \ 

\ X lM d 9 x Pe{x,a)4){x) H -E + /j,(e) J 
Then, one can construct the regularized resolvent of this hamiltonian using an explicit 
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formula given by Rajeev [T] 




where 



ME) = ^— -g = | (23) 



-Ho — -Ho — E Hq — E 

5 = *-\E) 

b = X d g x p e (x,a)(f)(x) . (24) 

Jm 

Here E should be considered as a complex variable. Most importantly, the operator & e (E), 
called principal operator, is given as 

§ € (E) = H - E + /i(e) - A 2 / d g x d g y p e (x, a)p £ (y, a) 4>{x) $(y) , (25) 

Jm ±i — h, 

Once we have a proper definition of the principal operator, all the divergences are removed 
since the resolvent is expressed in terms of it. We can extend the principal operator by 
analytic continuation to its largest domain of definition in the complex energy plane. For 
our purposes, we will assume that 31(E) < nm + p,. In fact, the energy of bound states 
are interesting and they satisfy the required conditions as we will see. Now, we will do 
the normal ordering of the operators in (|25p by using the commutation relations of the 
operators 4>(x) and (ft(x) 

0(x)—^— = [ d g x I'™ dse- s ^- E ^K s/2m { X , X i)4>{x'), 
no — n, j M j 

which can be proved simply by using an eigenfunction expansion. Then, the principal 
operator can be written in terms of heat kernel 



poo P 

$ e (E) = H -E-\ 2 ds d g xdgydgx'd g y'p e (x,a)p e (y,a) 

Jo Jm 

x K s/2m (x, x') K s/2m (y, y') </>V) e -< H " +2m -^ <t>(x 



poo 

+ /i(e)-A 2 / ds I dgyd g y f K s/2m (y,y')p e (y,a)p t (y',a)e~ s( - Ho+m - E K 
Jo Jm 

Since the heat kernel is a natural delta sequence, we can set p e (x, a) = K e u m (x, a) without 
loss of generality. Hence, 

poo p 

$ e (E) = H -E-X 2 ds dgxd g ydgx' d g y' K e/Am (x,a)K e/im (y,a) 

Jo Jm 

x K s/2m (x, x') K s/2m (y, y') <j>\y') e -s(H 0+ 2rn-E) ^ 



poo p 

+ /i(e) - A 2 / ds d g yd g y' K s/2m (y, y')K t/im {y, a)K e/4m {y\ a) e -s(H 0+m -E) 
JO J M 

Using the reproducing property of the heat kernel, one can rewrite the above equation as 

poo p 

$ e (E) = H -E-X 2 ds d g xd g y K(2s+e)/4m(x,a)K(2s+e)/4m(y,a) 

JO JM 

poo 

X <p(x) e s(H 0+ 2m-E) m + M(e) _ A 2 / ds K{s+e)/2m{a: a) e s(H 0+m -E) 



Shifting the variable 5 in the first integral by e/2 and in the second integral by e, we get 

pOO p 

$ e (E) = H -E-X 2 ds d g xd g y K s/2m (x,a)K s/2m (y,a) 

Je/2 JM 



oc 



X <fi{x) e -(s-*m(Ho+2rn-E) ^ + ^ _ y2 J ^ K s/2m (a, a) e -(s-e)(H 0+ m-E) _ 

If we take the limit e — > + , the only divergence is coming from the lower limit of the 
second integral term due to singular behavior of the diagonal part (k = term in the 
sum) of heat kernel near s = 0, see equation (|lip . 

One can also see that the first interaction term is actually finite due to the quite sharp 
bounds on the heat kernel for various classes of manifolds [HJ [13] . In fact we will explicitly 
show later on that this term is really a finite expression by working out a bound on the 
spectrum of the model. Since the principal operator or resolvent is not well defined in 
this limit, we must now regularize the model by choosing [i{e) exactly the same as in (|13|) . 
Then, we find 



x </>t 



/"CO p 

$ t (E) = H - E - X 2 ds d g x d g y K s/2m (x, a)K s/2m (y, a) 

Je/2 JM 

/oo 
dsK s/ 2 m (a,a) 



e -s{m-ij) _ e ~s(H a +rn-E) 



e -s{m-fi) _ e -s(H +m-E) 



Here the limit e — > + is now well-defined so we have 

poo 

= H -E + n + X 2 dsK s/2m (a,a) 
Jo 

poo p 

-A 2 / ds dgxdgyK s/ 2 m (x,a)K s/ 2 m (y,a) ( l ) \x)e~ s( - Ho+2m - E U(y) ■ (26) 
Jo Jm 

This is the renormalized form of the principal operator so that we have a well-defined 
explicit formula for the resolvent of the Hamiltonian in terms of the inverse of the principal 
operator $ _1 (i?). 

The spectrum of the Hamiltonian is the set of numbers E at which the resolvent does 
not exist (discrete spectrum) or exist but is unbounded (continuous spectrum). Thus, 
the continuous spectrum is that of Hq and the values of E where $(£J) does not have a 



bounded inverse. Since, there are no poles in H ^_ E , the poles corresponding to bound 
states must arise from those of ^^ 1 (E), that is, the roots of the equation 

$(£)|tf)=0, (27) 

determine the poles in the resolvent, which means that the principal operator &(E) de- 
termines the bound state spectrum of the theory. After we have found a root, we can 
determine the corresponding eigenstate of the Hamiltonian. A trivial example is the 
bosonic vacuum state, 

d>(£)|0> = jff - E + fj, + A 2 ds K s/2m (a, a) [ e -»(^) - e -s{H 0+m -E) 

-A 2 f°° ds [ d g xd g y K s/2m (x,a)K s/2m (y,a)^(x)e-^ Ho+2m - E U(y))\0) = , 



'0 JM 

where the root can be easily found to be 

E = Ll. 

We remark here that the linear eigenvalue problem is converted into a nonlinear problem 
for an operator family parametrized through the energy eigenvalues after the renormal- 
ization. 

The residue of the pole in the resolvent is the projection operator to the corresponding 
eigenspace of H [14] 

where is a small contour enclosing the isolated eigenvalue fx in the complex plane. For 
this contour integral we need to evaluate the residue 

Res^^-^IOXOI) = ($'(^))- 1 |0)(0| . (29) 

As a result of this calculation, we find the projection operator to be 

poo "1 —1 

— 1 i \2 / „ ^„ „-s(m-n) v . („ „\ 




(30) 



1 + A 2 / sdse-< m -riK s/2m (a,< 

0t(o)|O)(O| 

|o)(o|^(«)t^ |o><o| 

Thus, one can read off the eigenvector of H corresponding to the root E = \i from 
the projection operator and then find the normalizable eigenstate (|16p with the correct 
normalization factor f|18[) . This eigenstate (|16p is the first excited state (eigenstate of 
neutron), not the vacuum in the whole Hilbert space B<S>C 2 . Also, it is easy to see that 

in 



the zero eigenvalue of the Hamiltonian corresponds to the proton state I I . Here 

V / 

m > fx for bound states since we want this model to describe the attractive interaction of 
such a two state system with bosons. However, it is not clear whether the proton is the 
state of lowest energy. There may be states which contain many bosons that have a lower 
energy. These questions will be answered in studying the principal operator as we will do 
in the next section. 

One can also generalize these ideas into the relativistic regime with coupling constant 
renormalization, but this needs further investigations so we are not going to elaborate on 
this and will only consider the non-relativistic Lee model. 



3 A Lower Bound on the Ground State Energy 

After the renormalization of our model, we must look at the spectrum of the problem 
because there are many theories in which even after the renormalization there are still 
divergences that makes the spectrum not bounded from below pQ. In this section we will 
restrict E to the real axis. In order to give the proof that the energy E is bounded from 
below, following the same idea as in [I], we first split the principal operator as 

$(E) = K{E) - U(E) , (31) 

such that 

K(E) = H -E + fi, (32) 

and 

e -s(m-fi) _ e -s(H +m-E) 
pec ~ p 

+A 2 / ds d g xd g yK s/2m (x,a)K s/2m (y,a)^(x)e-^ Ho+2m - E U(y)-m 
JO J M 

It follows immediately that K(E) > nm — E + /i, so it is a positive definite operator from 
our assumption E < nm + fi. Due to the positive definiteness of heat kernel K s / 2m (a, a ) > 
0, the difference of the two exponentials is a positive operator. As a consequence of this, 
the first integral term in U(E) is a positive operator and we can claim that 

U(E) < U'{E) , 

where 

f'OO f' 

U'(E) = A 2 / ds d g x d g y K s/2m (x, a) K s/2m (y, a) ${x) e -s(H 0+ 2m-E) 
JO J M 
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U(E) = -A 2 / dsK s/2m (a,a) 



This clearly forces 

9(E) > K(E) - U'{E) , 

or rewriting it as 

9(E) > K(Ef/ 2 (l - U'(E)) K(E) 1 / 2 , (34) 

where U'(E) = K(E)~ 1 / 2 U'(E) K(E)~ 1 / 2 and K(E), U'(E) are positive operators (so 
is U'(E)). We will now show that by choosing E sufficiently large enough it is always 
possible to make the operator 9(E) strictly positive, hence it is invertible, and has no 
zeros beyond this particular value of E. Therefore, if we impose 

\\U'(E)\\ <1, (35) 

then the principal operator 9(E) becomes strictly positive. In order to do some estimates, 
we will rewrite the interaction term in terms of eigenfunctions f a (x) and shift the operator 
$(x) to the leftmost side and the operator (j)(x) to the rightmost side 

U'(E) = \ 2 Y1 ^Vi) UM Wo + /u + a x j2m + m - E]~ 1 ' 2 

0"1,0"2 

x [H + a 1 /2m + a 2 /2m + 2m - E]' 1 ^ + n + a 2 /2m + m - E}~ 1/2 4>(a 2 ) / CT2 (a) , 

where 4>(a) = f M d g x f*(x)4>(x). In order to convert the product of operators in the 
above formula into a summation of them, we will use the Feynman parametrization [15] 



1 1 L . A 1 „ . A 1 „ . A 1 rvi -4-rvo-l-rvQ ' V ' 



A^A^Af 1(^)1 (a 2 )T(a 3 ) Jo [«i^i + u 2 A 2 + u 3 A 3 }^+^+^ 

so that, 

fnim !iVA w i„\ £(1/2+1/2+1) 
u(E) = „^ r(i/2)r(i/2)r(i) 

«1 



x I dui du 2 du 3 u'l u| u 3 $( u l + u 2 + u 3 — 1 

1 







[flo + m + n - E + ( Ul + u 3 )^ + (u 2 + u 3 )§± + u 3 (m - fj,)} 2 
or one can rewrite it as 

T(2) f 1 du\ du 2 du 3 5(u\ + u 2 + u 3 — 1) 

0"1,0"2 



V2)fa 2 (a) , 



?/'(£) = A 2 ^0t (ai)/ai(a) J_^_^ 



(nm 2 )V2 

X / S ds e - s (- H "o+m+At-£;+U3(m-/i)+(«i+«3)((T 1 /2m)+(u2+-"3)(o'2/2m)) f U2 (a) 







Let us express this equation in terms of the heat kernel: 

A 2 T(2) f°° f f 1 dm du 2 du 3 5(u x + u 2 + u 3 - 1) 



^ - ww[ sds ! M daxd9y L 



( UlU2 y/ 2 



1? 



X <f>\x) ^ s (Mi+U3)/2m( X ! a ) ^s(u 2 +U3)/2m(2/i 



a) e 



-s(Ho+fj,+m-E) —su 3 (m—fi) 



It is easy to see that heat kernel has the following scaling property in 3 dimensions: 

K a 2 s (x, y; g) = a~ 3 K s (x, y; oT 2 g) , (37) 
where g — ► a~ 2 g means that the metric gij is scaled by a conformal factor a~ 2 . Thus we 



have 
U\E) 



A 2 T(2) [°° f 

r(i/2) 2 J s J M a ( u i+ u 'i)' 1 9 xa (u2+u 3 )- i gy 

11 du\ du2 du 3 5{u\ + U2 + u 3 — 1) 



<p ] (x)K s/2m (x,a;(u 1 + u 3 ) x g) 



K s/2m {y, a; (u 2 + u 3 ) 1 g) e 



-s(H +n+m-E) -su s (m-n) 



<Kv) 



In addition, under the scaling of metric, the commutation relations obey the following 
rule 

[<P a ig{x)A ] a 2 g {y)] = d a 2 g (x,y) , (38) 



or 



[a 3/2 a 2 5 (x),a 3/2 ^ 29 (y)] = 6 g (x,y) . 



(39) 



This lead us to find the scaling property of the creation and annihilation operators 
4> a 2 g (x) = a~ 3 / 2 <f> g (x). Using <p( Ul +u 3 )- 1 g( x ) = ( n i + u 3) figfe) we define the creation 
and annihilation operators with respect to the new metric and obtain 



U'{E) 



A 2 r(2) r 00 



r(i/2) 2 J 



sds I d 

M 



(u 1 +u 3 )~ 1 g x d( U2 +u 3 )- 1 gy 



1 du\ du2 du 3 5(ui + U2 + u 3 — 1) 



/o (uiu 2 ) 1 / 2 (u 1 + u 3 ) 3 / 4 (u 2 + n 3 ) 3 /4 
x K s /2m {x, a; (ui + u^g) K s/2m (y, a; (u 2 + u 3 )~ g) 

X iu 1+ u 3 )-^) e- s ^+^- £ ) e-^-»U(u 2+U3) -ig(y) ■ 

In order to give an upper bound estimate on the norm of the operator U'(E), we apply 
the Cauchy-Schwartz inequality in the norm corresponding to the new metric and get 

A 2 T(2) f°° e _ s{nm+fl _ E) f 1 du 1 du2du 3 5(u 1 +U2 + u 3 -l) 



\U'(E)\\ < n 



r(l/2) 2 



M 



d {ui+U3) -i g x K 2 /2m (x, a; (m + u 3 ) l g) 



(n 1 n 2 ) 1 /2( Ul+U3 )3/4 (u2 + ?/3) 3/4 
1/2 



M 



d( U2+U3 )-i g y K 2 /2m (y, a; (u 2 + u 3 ) 1 g) 



1/2 
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Here we have replaced the term Hq + m + \x — E in the exponent by its minimum value 
ran + fx — E and dropped the term e - su 3( m ~v) < \ gy us i n g the reproducing property of 
the heat kernel, we get 



< "Hr-'"-'"'-'!' 



A 2 T(2) f°° nJ __ s(nm+fl _ E) f 1 d Ul du 2 du 3 

:<5(«i + u 2 + n 3 - 1) [A' s/m (a, a; (ui + u^g)] 1/2 [^ s / m (a, «! («2 + ^3) _1 ff)] 1/2 • 



Scaling back to the original variables, we finally obtain 

x [K s(ui+U3)/m (a,a;g)] 1/2 [K s(u2+U3)/m (a,a; g)] 1/2 . (40) 

It is essential here to note the behavior of the heat kernel. In most of the situations, 
the explicit form of the heat kernel is unknown so one should look for the estimates or 
upper bounds on it. Luckily, there are quite sharp upper bounds on the heat kernel for 
various classes of manifolds in the mathematical literature [8l [T3"l [T6l [T8] . For each class 
of manifolds, there are different bounds so we will consider them separately. 

We will first consider Cartan-Hadamard manifolds, which are geodesically complete 
simply connected non-compact Riemannian manifolds with non-positive sectional curva- 
ture — K 2 (for example W 1 and M. d ). On 3 dimensional Cartan-Hadamard manifolds, we 
have the following upper bound on the heat kernel (see chapter 7.4 in [8] and [16j ) 



K s/2m (x, x) < . 1 ° e -*V»n , (41) 

or for simplicity of our calculations, one can use the following less sharp estimate [8] 

K s/2m (x,x) ,< (42) 



for all x G Ai and s > 0, and C is a dimensionless constant. This estimate (|42h can be 
extended also to the minimal submanifolds of M 3 [8]. Then the strict positivity of the 
principal operator after taking the integrals leads to the following inequality 

\\U'(E)\\ < n A 2 C H + M - £)- 1/2 < 1 • 

This implies the inequality for E or the opposite inequality for the ground state energy 
written in ordinary units 

2 2 2 C 2 A 4 m 3 „ 2 2 oC^AV 3 

£ < nmc + /ic — n — ^ — , E gr > nmc + Lie — n — ^ — , (43) 
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where 

- vrV 2 r(2)r(l/4) 2 

r(i/2) 2 r(3/4) 2 ' 

Secondly, as an explicit application, we consider the 3 dimensional hyperbolic space H 3 . 
The heat kernel in H 3 is explicitly known [U] and it is in natural units 



. . _ 1 d(x,y)/R e 2ms z a s 

where R is the length scale and d(x, y) is the geodesic distance between two points on H 3 . 
Again we impose the strict positivity condition on the principal operator and get 

A 2 T(2) f 1 d Ul du 2 



T(l/2) 2 J ( Ul n 2 )V2(i _ Ul )3/4 (1 _ U2 )3/4 

poo ^ s e -s(nm+fi-E) e -s(l-ui)/2mR 2 e -s(l-u 2 )/2mR 2 

Using the fact that e~ s ^ l ~ Ull - 2 ^^ 2mR < 1 and then taking the integrals, one can immedi- 
ately see that this implies 

^ ^ 2 2 2 A 4 D 2 m 3 ,. e , 
hi gr > nmc + [j,c — n -g , (45j 

where the constant D is defined as 

r(2)r(i/4) 2 



D 



r(i/2) 2 r(3/4) 2 2 3 ' 

Finally, we apply our method to the closed compact manifolds with Ricci curvature 
bounded from below by — k. The estimate for this class of manifolds on the heat kernel 
[T3"l [18] is given by 

K s/2m (a, a) < + A(s/2my 3 / 2 , (46) 

where A = A(V(Ai), n, d) is an explicitly calculable constant which depends on the volume 
V(M) of the manifold, the lower bound n on the Ricci curvature and the diameter d of 
the manifold. Using this estimate, it follows that 

\\U'(E)\\ < n-^t r 8dae -«™+^*) f 1 du x du 2 



r(i/2) 2 7 J {mm) 1 ' 2 



V(M) 1 / 2 



Integrating with respect to u±, u 2 and s, we have 



\U'(E)\\ < n 



X 2 T(2) 
r(l/2) 2 



4 1 4^ 1 /2 m 3/4 7r l/2 r(1/4) 



+ 



V (M) (nm + [i - E) 2 V l l 2 T{?>/4) 
IF) 



1 Am 3 / 2 vrr(l/4) 2 1 



(nm + /i-£) 5 / 4 T(3/4) 2 (nm + n - E) 1 / 2 _ 

In order to get explicit estimates, let us put some further assumption nm + n — E > fi. 
Then, we find 



\U'(E)\\ < n 



A 2 r(2) 



4 4A 1 /2 m 3/4^1/2 r ( 1 / 4 ) Am 3 / 2 7lT(l/4) 2 



r(i/2) 2 

1 

X 



+ ^77Z ' ' + 



V(M)a s / 2 f i^V(M) 1 / 2 T(3/A) r(3/4) 



(nm + u — E) l l 2 

Now if we impose the strict positivity of the principal operator we find similarly 

E gr > nmc 2 + ac 2 - n 2 X 4 F 2 , (47) 



where 



F 



r(2) 
r(i/2) 2 



4^1 1 /2 m 3/4 7r i/2 r ( 1 / 4 ) Am 3 / 2 7rr(l/4) 



+ ./ + 



y(M)/x 3 / 2 fi 3 / 4 V(My/ 2 T(3/4) T(3/4) 2 
Therefore, lower bound on the ground state energies for different classof manifolds (1431) . 
(|45p and (|47p are of almost the same form up to a constant factor so the form of the lower 
bound has a general character. It is also worth pointing out that the form of the lower 
bound on the ground state energy is same as in the case for the flat space M 3 [H I17j . 
From the general form of the lower bounds, we conclude that for each sector with a fixed 
number of bosons, there exists a ground state. However, the ground state energy bound 
that we have found diverges quadratically as the number of bosons increases. In other 
words, these estimates with our present analysis is not good enough to prove the existence 
of the thermodynamic limit. To attack this problem we will study the thermodynamic 
limit of the model by mean field approximation. 



4 Mean Field Approximation of the Model 

In the limit of large number of bosons n — > oo, one expects that all the bosons have the 
same wave function u(x) and mean field approximation is valid, as in the case of flat 
spaces. Therefore, one can introduce the following mean field ansatz for n particle state 
on a Riemannian manifold 

\u) = — f= / dgxi d g x 2 ■ ■ -dgx n u{x\) ufa) ■ ■ ■ u(x n ) 4>\xi) (j)\x2) ■ ■ ■ ^(x n )\0) , (48) 
Vnl Jm 

with the normalization 

\\u(x)\\ 2 = / |n(x)| 2 d g x = 1 , (49) 

Jm 

if? 



where |0) denotes the vacuum state. In the mean field approximation the operators are 
usually approximately replaced by their expectation values in this state ((/(«)) ~ /((«))) 
so the expectation value of the principal operator by applying the mean field ansatz 
becomes 

/•CO 

4>(E, u) = nh (u) -E + fi + X 2 ds K s/2m (a, a) [ e -< m ~^ - e -s(nh (u)+ m -E^ 

Jo 

/•OO P 

- nX 2 / ds / d g x d g y K s/2m (x, a) K s/2m (y, a) u*(x) e -<nh {u) + 2 m -E) u{y) j (5Q) 
Jo J M 

called principal function and we have defined 

h (u)= [ d g x ( IV ^ {X)]2 +m\u(x)A . (51) 

Jm V 2m ) 

However, the exact value of the expectation value of the principal operator is given in 
terms of cummulant expansion theorem [19]. Therefore, in order to write the above 
formula (|5U|) . we have to further assume that the corrections coming from the higher 
order cummulants are negligibly small and indeed we will see that this assumption is 
justified for the particular solution we will find. 

Now, we must solve the equation 4>(E, u) = (due to (|27p ) in order to find the spectrum 
of the problem and solve £ as a function of u(x), which gives the smallest possible value 
of E with the constraint (|49p . Hence, one can try to write E as a functional of u(x) from 
the equation (p(E, u) = and apply the variational methods to minimize E. 

However, there is no simple way to solve this variational problem . So, we follow a 
different method essentially the one suggested in pQ. For this purpose, let us introduce a 
new variable x an d new wave function v(x) 

X = nh (u) - E 
v(x) = [2m(2m + x)]~ 3/4 u(x) , (52) 

such that the new wave functions v(x) are normalized with respect to the new metric 
gij = [2m(2m + x)]9ij 

I d g x \v(x)\ 2 = I d g x \u(x)\ 2 = 1 . (53) 
Jm Jm 

We also define a new dimensionless parameter s' = (2m + x) s i an d using the scaling 
property of heat kernel we find 



X + n + A 2 / ds K s/2m (a, a) 
Jo 



e -s(m-fi) _ e -s(x+m) 
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f'OC r 

nX 2 (2mf /2 (2m + x) 1/2 / ds 1 / d g x K s ,{x, a; g) v( 

Jo Jm 



(r) 



(54) 



One can prove now that the left hand side as a function of the variable x is an increasing 
function and it is obvious that the righthand side is positive. Therefore, the left hand side 
is minimum when x = which is attained when the right hand side becomes zero (so 
X > — A*)- Let us denote the inverse function of the left hand side as fi(nU), that is, 

X = fi(nU), (55) 

and express x i n terms of the energy E and the function v (x), 

X = n[x + 2m]i^[?;] + nm — E, or (56) 

E = nm + 2mnK[v] + (nK[v) - l)fx(nU) , (57) 

where K[v] = f d g x\Vgv(x)\ 2 , which is considered to be the parameter of the model 
because it is the variable we can control and we may use many trial functions v(x) and 
they can be scaled to any desired value. If we assume that ni^[w] > 1, then the energy 
E (|57p is minimized when fi(nU) is minimized which happens for U[v] = 0. Since 
X = fi(nU) > — /i, we have 

E>nm + fi if nK[v) > 1 . (58) 

On the other hand, if K[v] is small enough, i.e., raET[u] < 1, we also see that the minimum 
of the energy is attained with the reversed sign of the last term. In that case, we should 
find an upper bound for f\{nU) which is expressed in terms of the kinetic energy functional 
K[v]. In order to discuss the case n-fffv] < 1 properly, we will separate our calculations 
for compact and non-compact manifolds. 

Let us first consider the case for compact manifolds. In order to achieve our aim, we 
will go back to the original variable u(x) and the parameter s in the equation (154j) 



roc 

X + + A 2 / ds K s/2m (a, a) 
Jo 



e -s(m-fi) _ e -s{x+m) 



n\ 2 I ds 







d g x K s/2m (x,a;g)u(x) 
M 



2 



e 



If we assume that the mean field approximation gives us a reliable equality, we can find 
the solution for x given all the other parameters. Note that the right hand side of (|59p is 
a decreasing function of x whereas the left hand side is an increasing one. Hence there is 
always a unique solution which defines the inverse function for a given u(x). It is quite 
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easy to see by a graphical construction that if we replace the left hand side by a smaller 
function of x> an d the right hand side by a larger function of x> then we get an upper 
bound for the inverse function. We note first that 

/ d g x K s/2m (x,a;g)u(x) = V) e~^/ CT (a) u(a) , 

Jm 

where u(a) = J M d g x f£(x) u(x). If we define the following functions 



( 2mf a (a) 



/» = 



and 



{ Ha) 
au(a) 



u\a) = < 



2m 



u(0) 



if a^O 



if a = 0, 



if a ^ 



if a = 0, 



(60) 



(61) 



we can write 



M 



d g x K s/2m (x,a;g)u(x] 



2m /£(a)u' \a) 



<E e ""l/>)! 2 El-V)| 2 < (l/o(a)| 2 + 2m^^|/ CT (a)| 2 j ( |u(0)| 2 + 



K[u] 
2m 



where is the sum which excludes the zero mode (a = 0), K[u] = ^ j M d g x\ V g u(x)\ 2 , 
and we have used Cauchy-Schwartz inequality. Using |/o(a)| 2 = 1/V{M) and |w(0)| 2 < 1, 
we find 

u[u] < (i + n , (62) 

where 

A 2 



n 



X + 2m 



i r°° 



x [K s / 2m (a,a;g) - lim K s/2m (a, a; g) ) 



• (63) 



Here the sum which excludes the zero mode corresponds to substraction of the large s 
behavior of heat kernel. Using K[u] = (2m + the inequality becomes 



X + M < 



X + 2m 



x ( K s/2m(a,a;g) - lim K s/2m (a, a; g) 



i + (x + 2m) ^] 

2m L J 



If) 



Using the upper bound estimate of the heat kernel for closed compact manifolds with 
Ricci curvature bounded from below by —k and taking the integral with respect to s, we 
find that 



X + M < 



nX 



2 



X + 2m 



1 



V{M) 



+ V2^{2mf/ 2 A( X + 2m) 1 / 2 \ (64) 



If we now introduce the variables z = \ + 2m and A' = y/2ir(2m) 5 / 2 A for simplicity of 
notation, we find the following inequality 

,_ (2m _,)<^(„ + -i-)(^ + _±_). (65) 

We look for a systematic expansion of z in n by allowing a fractional power. In this 
case we see that if we substitute the following asymptotic expansion as n — ► oo, z ~ 
B\n Vx +B2n u ' 2 + B^n V3 + • • • , where the consecutive powers u±, . . . decrease, we find that 
this asymptotic expansion is an upper bound for this inequality as long as the coefficients 
are to be chosen as B l = {A'X 2 ) 2 / 3 , B 2 = (jpf /3 + ^(A'X 2 )^ 3 and B 3 = 2m - 

/i + 2mV{M) ' wnere ^1 = 2/3, ^2 = 1/3, and z/3 = 0. Hence we get an upper bound on the 
inverse function f\{nU) as an asymptotic series in powers of n, in the spirit of mean field 
approximation. This upper bound can be put back into the energy equation and then we 
find the following lower bound on the ground state energy, 

E gr >nm + 2mnK[v] - (1 - nK[v\) {Bin 2/Z + B 2 n 1/3 + B 3 - 2m H ) (66) 

We note now that the behavior of K[u] for large n is found from the scaling law, K[u] = 
(2m + ~ ra^ 1 / 3 . This in turn justifies our use of the mean field approximation 

since higher order derivatives are then negligible, and the solution is approaching to an 
essentially constant function on the manifold as n gets larger. This proves that for a 
compact manifold the energy is actually bounded form below by a much milder behavior 
and there is a nice thermodynamic limit since the energy per particle E/n will approach 
to the mass (rest mass energy) as n — ► 00. The same conclusion can also be drawn for 
non-compact manifolds as we will see. 

Now we will consider the mean field approximation of the model for non-compact 
manifolds. We again assume that the eigenfunction expansion in compact manifolds can 
be generalized to the non-compact manifolds. Then, one can use the above method for 
the non-compact manifolds as well. However, we shall try to find the ground state energy 
in the mean field approximation in an another way. Therefore, we first go back to the 



9,n 



equation: 



\2 



( X + 2m)^\ X + ^ + X Jo dsK °^ a) 



e -s(m-/i) _ e -s(x+m) 



nX\2mf/ 2 / ds' 
Jo 



M 



dgx K s ,(x,a;g)v(x) 



n U[v] , 



and using the generalized eigenfunction expansions, we find 



u[v) < k[v] n 



(67) 



where 



n 



A 2 (2m) 3 / 2 ds' (l - e^ s ' /2 ) K s ,(a, a; g) , 



(68) 



The explicit form of the inverse function f\{nU) is too difficult to find so one can estimate 
it. In order to do this, let us first notice that 

/ 2 _1 (X) < fiHx) then h(nU) < f 2 (nU). 

For this purpose, we use the simplest possible function as /^ 1 (x) 

f -i, \ X + V 

12 {X> (x + 2m)V2 • 

Then we can replace f\{nU) with something bigger, and its argument with something 
even bigger. Moreover, f2{u) is dominated by a simpler function 

f 2 (u) < u 2 + 2m - 2fi, 

which is a very crude bound, but easy to work with. Using the upper bound for U[v], we 
get 

x = ft(nU) < f 2 (u) < n 2 U[v] 2 + 2m-2fi< n 2 (K[v]) 2 n 2 + 2m-2^. (69) 

Hence, 

E>nm + 2mnK[v] - (1 - nK[v\) {n 2 K[vf Q 2 + 2m - 2/i) , (70) 

where the polynomial (4m — 2\i)y — Q 2 y 2 + 2 y 3 never becomes negative within the range 
< y = nK[v] < 1 if 

n 2 < 16m - 8/i . (71) 

In this case, the minimum is achieved when y = 0. This means that the functional 
nK[v] — > by a proper family of functions. This is why it is consistent to ignore the 
higher order cummulants if we choose an arbitrarily slowly varying family for the function 
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v{x). Therefore one can conclude that the ground state energy for non-compact manifolds 
when nlfkv] < 1 

E gr > nm — 2(m — /i) . (72) 

The condition on O may imply some restrictions on the coupling constant A. If we go 
back to the definition of and scaling back again to the usual geometric variables, we 
find 

n = \ 2 ( X + 2m)- 1 / 2 J™ds (l-e-<* +2m V 2 )K s/2m {a,a;g) . (73) 

Since there is a nice sharp upper bound on the heat kernel for Cartan-Hadamard manifolds 
and minimal submanifolds of K 3 (|42p . we can find an upper bound on 17 and the restriction 
(I71|) gives the following upper bound on the coupling constant. 

A < {16m ~ 8/i)V4 (74) 
(27r) 1 /4cV2(2m)3/ 4 ■ 1 ' 

Now, let us calculate explicitly the function Q for the hyperbolic manifold H 3 , which 
belongs to the class of Cartan-Hadamard manifolds. Using the result (|44p . we obtain 



Then, one can easily find the upper bound on the coupling constant from the restriction 
on 0, 

A < 2VM2m - ,){2my^{j 2 -^ + - JZ- .CO \ (76) 



Therefore, similar to the case for compact manifolds, we have shown that the leading 
behavior of the system varies linearly with the number of bosons n, which leads to a nice 
thermodynamic limit on non-compact manifolds. 



5 Conclusion 

In this paper, we considered the non-relativistic Lee model on various class of Riemannian 
manifolds inspired from the work in pQ. This method allows us to renormalize the model 
non-perturbatively. It has been also shown that the heat kernel plays a key role in the 
renormalization procedure and help us to find a lower bound on the ground state energy 
due to the sharp bound estimates on it for several class of manifolds. Finally, we studied 
the mean field approximation and showed that there exist a nice thermodynamic limit of 
the model for compact and non-compact manifolds. 



9,9. 
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